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Abstract. Both high resolution spectra of QSOs observed at the 8-10 m 
telescopes and advanced methods of data analysis are crucial for accu- 
(j ' rate measurements of the chemical composition and physical parameters 

. of the intervening clouds. An overview of our recent results obtained with 

the Monte Carlo inversion (MCI) procedure is presented. This includes: 
(1) variations of the shape of the local background ionizing continuum in 
the 1-5 Ryd range at redshift z ~ 2.8 — 3.0; (2) an inverse correlation 
between the measured metallicity, [C/H], and the absorber line-of-sight 
linear size, L; (3) a functional dependence between the line-of-sight ve- 
^ ■ locity dispersion, a v , and L. 
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1. Introduction 

A study of quasar (QSO) absorption-line spectra is generally recognized as the 
most reliable technique for inferring the physical and dynamical state of gas in 



the intervening absorption clouds at high redshifts, z ;> 2. Of particular interest 
are the measurements of metallic absorptions in the optically thin diffuse clouds 
with neutral hydrogen column densities N(R i)^3x 10 17 cm -2 . These systems 
will be called as 'Ly-a absorbers' (LAA) to distinguish from the damped Ly-a 
absorption systems (DLA) which show much higher column densities, ^(H i) ^ 
2 x 10 20 cm~ 2 . The latter are believed to arise in the galactic disks (e.g., Wolfe et 
al. 1995), whereas the former may be physically related to the external (~ 10-100 
kpc-scale) regions of galaxies (e.g., Chen et al. 2001). Thus the measurements 
of the LAAs can provide fundamental insights into conditions prevailing in the 
galactic environments (external halos) in the early universe. 

Since these regions are mainly photoionized by the local metagalactic UV 
radiation, the ionization states of the LAAs are sensitive to the spectral shape of 
the background radiation in the 1-5 Ryd range. The spectral energy distribution 
in the metagalactic ionizing background is defined in turn by the QSO continua 
filtered through the quasar environments and the IGM. Two implications of the 
LAA analysis - the metal content of the external halos and the spectral shape 
of the local UV background - are, therefore, tightly coupled. 

To clarify the mechanism of the metal enrichment, accurate measurements 
of the metal abundances in the LAAs are required. The main problem here is 
how to account for the ionization correction. In general, the contribution to the 
line intensity I\ within the profile comes from all volume elements distributed 
along the line of sight and having the same radial velocity. If the gas number 



1 



2 



S. A. Levshakov 



density, bh, varies from point to point, then the intensity I\ is caused by a 
superposition of different ionization states. 

Recently, we have developed a method called 'Monte Carlo inversion' (MCI) 
to recover the physical parameters of the LAAs assuming that the absorb- 
ing cloud is a continuous region with fluctuating density and velocity fields 
(Levshakov et al. 2000, hereafter LAK). The MCI was applied to high qual- 
ity QSO spectra obtained with the VLT/UVES, Keck/HIRES, and HST/STIS 
(Levshakov et al. 2002; Levshakov et al. 2003a,b,c,d,e). The results of these 
studies are briefly reviewed in this contribution. 

2. The MCI procedure 

The layout of the MCI procedure is the following. We assume that the metal 
abundances within the absorber are constant, the gas is optically thin for the 
ionizing UV radiation, and the gas is in thermal and ionization equilibrium. 
The radial velocity v{x) and the total hydrogen density «h(s) along the line 
of sight are considered as two random fields which are represented by their 
sampled values at equally spaced intervals Ax, i.e. by the vectors {v\, . . . ,Vk} 
and {ni, . . . with k large enough (~ 150 — 200) to describe the narrowest 
components of the complex spectral lines. The radial velocity is assumed to 
be normal distributed with the dispersion a v , whereas the gas density is log- 
normal distributed with the mean uq and the dispersion a y {y = nn/no). The 
stochastic fields are approximated by the Markovian processes. The accuracy 
of the restoring procedure depends on the number of different ions in different 
ionization stages involved in the analysis of a given LAA. 

A set of the fitting parameters in the least-squares minimization of the 
objective function (Eqs. [29] and [30] in LAK) includes a v and a y along with the 
total hydrogen column density iVn, the mean ionization parameter Uo, and the 
metal abundances Z a for a elements observed in the LAA. With these parameters 
we can further calculate the mean gas number density no, the column densities 
for different species N a , the mean kinetic temperature T^, and the line-of-sight 
size L = N-ft/riQ of the absorber (note that uq and, correspondingly, L scales 
with the intensity of the local background radiation field) . 

We start calculations assuming some standard ionizing spectrum [e.g., power 
law, Mathews & Ferland (1987), or Haardt & Madau (1996, hereafter HM)] and 
compute the fractional ionizations and the kinetic temperatures at each point 
x along the sightline with the photoionization code CLOUDY (Ferland 1997). 
To optimize the patterns of {vi} and {rii} and to estimate simultaneously the 
fitting parameters, the simulated annealing algorithm with Tsallis acceptance 
rule and an adaptive annealing temperature choice are used. If the fitting with 
the standard spectrum is impossible, its shape is adjusted using the procedure 
based on the experimental design technique (Levshakov et al. 2003e). 

3. The spectral shape of the ionizing continuum 

In practice, the shape of the ionizing spectrum can be estimated in cases when 
the value of iV(H i) is measured accurately and the absorption system contains 
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Figure 1. A typical metagalactic ionizing spectrum at redshift z ~ 3 
(the dotted curve) modeled by Haardt & Madau (1996), and its modifi- 
cation (the solid curve) required to match the absorption lines observed 
in a LAA at z = 2.82 toward HE 0940-1050 (Levshakov et al. 2003e). 
The spectrum is normalized so that J v {hv = 1 Ryd) = 1. The emission 
bump at 3 Ryd is caused by reemission of He n Lya and two-photon 
continuum emission from intergalactic clouds. 



unsaturated lines of at least C Il-C iv and Si n-Si iv, otherwise the system is 
not sensitive to the fine tuning of the continuum shape. 

We estimated the spectral shape for the LAAs at z = 2.82 (H I, C II, C III, 
Si in, C iv, Si iv) toward HE 0940-1050; z = 2.7711 (H I, C II, Si II, C III, N III, 
Si in, C iv, Si iv, O vi), and z = 2.94 (H I, C n, Si n, C in, Si in, C iv, Si iv) 
toward Q 1157+3143 (Levshakov et al. 2003e). 

All three recovered spectra of ionizing radiation show common features: a 
bump at E = 3 Ryd, which is more pronounced comparing to the model mean 
intergalactic UV spectra at z = 3 like that of Haardt & Madau (1996), and a 
sharp break just after the bump - also at variance with the model predicting 
a smeared out break at E = 4 Ryd due to ionization of He II. An example 
of ionizing background estimated for the z = 2.82 absorption system toward 
HE 0940-1050 is shown in Fig. 1. This spectral shape rules out a considerable 
galactic contribution to the QSO dominated UV ionizing background at z ~ 3. 
The recovered UV spectrum can be well explained in the scenario of the delayed 
re- ionization of He n (Reimers et al. 1997). In this case the sharp break at 
E = 3 Ryd occurs due to strong resonant scattering of QSO radiation in metal 
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Figure 2. Carbon abundances [C/H] plotted against the logarithmic 
linear size L of the absorber estimated by the MCI procedure (Lev- 
shakov et al. 2002; 2003a,b,c,e). [C/H] decreases with increasing L 
reflecting, probably, the dilution of metals within galactic halos caused 
by the mass transport processes. 



and He n Lyman series lines whereas a part of the absorbed photons re-emitted 
by the intergalactic gas in the He n Lya and two-photon continuum emission 
increases the amplitude of the bump at 3 Ryd. Our results also indicate that 
the re-ionization of He n has not been yet completed by z = 2.77. This results 
is in line with recent observations of the He n Ly-a forest by Shull et al. (2003) 
who claim that 'the ionizing background is highly variable throughout the IGM' 
at z ~ 2.8. 



4. '[C/H]-L' and <a v - N n V relations 

The analyzed LAAs show that they are a heterogeneous population that is formed 
by at least three groups of absorbers : (1) extended metal-poor (Z < 0.1 Zq) gas 
halos of distant galaxies; (2) gas in dwarf galaxies (O.1Z < Z <^ 0.3 Zq); and 
(3) metal-enriched gas (Z ^ 0.5 Zq) arising from the inner galactic regions and 
condensing into the clouds within the hot galactic halo (high redshift analogs to 
the Galactic high velocity clouds, HVC). 
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Figure 3. Plot of the line of sight velocity dispersion log(<r„) vs. 
log(A r H L) for the same sample of the LAAs shown in Fig. 2. The 
dashed line corresponds to the relation \og(a v ) oc 0.5 log(A^H L) ex- 
pected for the virialized systems. Open squares represent HVC-like 
clouds with L ^ 1 kpc. 



Figure 2 shows a plot of the measured carbon abundances [C/H] 1 versus 
logarithmic sizes of the studied systems. Systematically higher abundances are 
seen in compact systems. This tendency reflects, probably, the dilution of met- 
als within galactic halos caused by the mass transport processes (like diffu- 
sion, turbulent mixing, galactic rotation, shear, and etc.). Metals are, probably, 
transported into the halo in form of small dense clouds carried out by wind or 
jets. In some cases we directly observe such blobs. For instance, the LAAs at 
z = 1.385 ([C/H] ~ -0.3, L ~ 1.7 - 2.5 kpc) and at z = 1.667 ([C/H] ~ -0.5, 
L ~ 1 kpc) toward HE 0515-4414 as well as that at z = 2.966 ([C/H] ~ -0.4, 
L ~ 100 pc) toward Q 0347-3819 are embedded in extremely metal-poor halos 
with [C/H] < -2 (Levshakov et al. 2003b,c). 

If LAAs are formed in gas clouds gravitationally bound with intervening 
galaxies, their internal kinematics should be closely related to the total masses 
of the host galaxies. We find a correlation between the absorber's linear size 
L and its line-of-sight velocity dispersion a v . The virial theorem states : a% oc 
M/L oc no L 2 = iVn L. Assuming that the gas systems are in quasi-equilibrium, 



1 Using the customary definition [X/H] = log(X/H) - log(X/H)o. Photospheric solar abun- 
dances are taken from Holweger (2001). 
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one can expect o v oc (N^L) 1 / 2 . In Fig. 3 we plot the measured values of o v 
versus the product of L and the total gas column density A%. It is seen that 
most systems with linear sizes L > 1 kpc lie along the line with the slope 0.5. 
Taking into account that we know neither the impact parameters nor the halo 
density distributions, this result can be considered as a quite good fit to the 
expected relation for the virialized systems. Hence we may conclude that most 
absorbers with L > 1 kpc are gravitationally bound with systems that appear to 
be in virial equilibrium at the cosmic time when the corresponding LAAs were 
formed. 
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